The structure of tin isotopes with a global optimized effective interaction 
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We present a systematic configuration-interaction shell model calculation on the structure of light 
tin isotopes with a new global optimized effective interaction. The starting point of the calculation 
is the realistic CD-Bonn nucleon-nucleon potential. The unknown single-particle energies of the 
ld 3 / 2 , 2si/2 and 0/in/2 orbitals and the T = 1 monopole interactions are determined by fitting to 
the binding energies of 157 low- lying yrast states in 102_132 Sn. We apply the Hamiltonian to analyze 
the origin of the spin inversion between 101 Sn and 103 Sn that was observed recently and to explore 
the possible contribution from interaction terms beyond the normal pairing. 
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I. INTRODUCTION 

Substantial experimental and theoretical efforts have 
been devoted in the past decade to explore the struc- 
ture of light tin isotopes [ll422| . The excitation ener- 
gies of the first 2+ states in Sn isotopes between 102 Sn 
and 130 Sn are established to pos sess an almost constant 
value (see, e.g., Refs. [HH Hi]). This may be under- 
stood from the simple perspectives of generalized senior- 
ity scheme [H, [H, [2(| and pairing correlation [13, [Hf ■ A 
more realistic description of these nuclei requires a knowl- 
edge of the effective interaction between the valence nu- 
cleons that govern the dynamics 0, [2!| . The complete 
wave functions thus calculated show large overlap with 
those of the generalized seniority scheme, especially for 
the low-lying states of isotopes close to the N = 50 and 
82 shell closures [H, [28| . The J = pairing channel of 
the effective interaction has been shown to play an essen- 
tial role in reproducing the spacings between the ground 
states and 2^~ states [28|. Possible deviations from the 
generalized seniority scheme were suggested from B(E2) 
measurements in the 2^ states Q. 

A microscopic shell-model description of the configura- 
tions of nuclei in the trans-tin region is a challenging task 
due to the scarceness of available experimental data and 
the near-degeneracy in energy of the relevant 0*77/2 and 
1^5/2 single-particle orbits 
model calculations of Refs. 
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In the earlier shcll- 
the spacing between 
the two orbits were taken to be £{0*77/2) — £(1^5/2) = 
0.2 MeV and 0.5 MeV. In Ref. Hi), excited states in 
103 Sn have been observed using in-beam spectroscopic 
methods. The measured spectrum of 103 Sn is very sim- 
ilar to those of 105 ! 107 ' 109 S n! with the spin-parity 5/2 + 
and 7/2 + for the ground state and first excited state. 
By a shell-model fitting procedure, the spacing between 
0*77/2 and ld 5 /2 single-particle orbits was predicted to be 
£(097/2) - £(1^5/2) = 0.11 MeV [H|. It has long been 
expected that the ground state spin of 101 Sn should be 
identical to those of 103 - 109 Sn 0] which can be approxi- 



mately viewed as one-quasiparticle states [30|. However, 
in Ref. [Hj], the configurations of the ground state and 
first excited state in 101 Sn were determined to be 0<7 7 / 2 
and ld 5 /2, respectively. The spins of these states are 
reversed with respect to those in 103 Sn. 



Shell model calculations with empirical interactions 
have been very successful in explaining the structure and 
decay properties of ligh t nuclei between 4 He and 100 Sn 
(see, e.g.. Refs. 34-t38j) and heavier nuclei around shell 
1361 . [39l 4C I . The key to these calculations is a 
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proper description of the monopole channel of the effec- 
tive interaction |4lj . which determines the bulk proper- 
ties of the effective interaction and governs the evolution 
of the effective single-particle energies (the mean field) 
function of valence neutron and proton numbers. 
The contribution of the monopole interaction becomes 
much more important with increasing valence nucleon 
numbers N since it is proportional to N(N — l)/2. The 
light tin isotopes between shell closures N = 50 and 82 
are the longest chain that can be reached by contempo- 
rary shell model calculations. They may provide an ideal 
test ground to study the competition between different 
terms of the monopole interaction. 



Realistic effective interactions obtained from free 
nucleon-nucleon potentials provide a microscopic founda- 
tion to shell model calculations (29|. Extensive previous 
shell-model calculations tend to suggest that the realis- 
tic interaction can give a satisfactory description of the 
multipole part but not the monopole channel, which may 
be due to the lack of three-body forces. The monopole 
interaction has to be treated empirically [421 ]. Thus we 
are motivated to fine-tune the monopole part of the real- 
istic interaction by fitting to available experimental data 
in tin isotopes. One may also get a limitation on the un- 
known single-particle energies of the orbitals le?3/2, 2si/2 
and 0/i n / 2 . We expect that the refined effective Hamil- 
tonian will give a better understanding of the structure 
of trans-tin nuclei. 
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II. MODEL SPACE AND OPTIMIZATION OF 
THE EFFECTIVE HAMILTONIAN 

We assume the doubly- magic 100 Sn as the inert core. 
For the model space we choose the neutron and proton 
orbitals between the shell closures N = Z = 50 and 82, 
comprising 0g 7 / 2 , W5/2, 1^3/2, 2s 1/2 and 0/i u / 2 . Wc 
also assume isospin symmetry in the effective Hamilto- 
nian. A common practice in full configuration interac- 
tion shell model calculations is to express the effective 
Hamiltonians in terms of single-particle energies and two- 
body matrix elements numerically (see, e.g., the Oxbash 
Hamiltonian package [HI), 

H = ^e a N a 

a. 

+4 zZ ^^p\ V \h3s)jTA ] j T . jajfj Aj T -j 5 ^,{l) 

afiS-fJT 

where a = {nljt} denote the single-particle orbitals and 
e a stand for the corresponding single-particle energies. 



ator. (j a jp\V\j-yjs}jT are the two-body matrix elements 
coupled to good spin J and isospin T. Ajt (^4jt) * s the 
fcrmion pair annihilation (creation) operator. For the 
model space we have chosen the effective Hamiltonian is 
such that it contains five single-particle energies and 327 
two-body matrix elements. Among the two-body ma- 
trix elements there are 167 elements with isospin T = 
and 160 elements with T = 1. The number of matrix 
elements of a given set of J and T is given in Table 
UJ The monopole interaction is defined as the angular- 
momentum-weighted average value of the diagonal ma- 
trix elements (jajp\V\j a jp)jT for a given set of j a , jp 
and T [4l|, |42[. For the chosen model space there are 15 
T = (and 1) monopole terms. 

The single-particle energies are assumed to be the same 
for all nuclei within the model space. They are given rel- 
ative to the neutron Of/7/2 state. The energy of the 1^5/2 
is taken as e(ld 5/2 ) = 0.172 McV [13]. The energies of 
other states have not been measured yet. They are ad- 
justed to fit the experimental binding energies of tin iso- 
topes. The single-particle energies of the proton orbitals 
are assumed to be the same as those of the neutron. 

The starting point of our calculation is the realistic 
CD-Bonn nucleon-nucleon potential [Hj|. The interac- 
tion was renormalized using the perturbative G-matrix 
approach, thus taking into account the core-polarization 
effects [29( . This interaction has been intensively applied 
in recent studies 0, Hf| 511 ■ The mass dependence of 
the effective interaction is not considered in the present 
work. 

Calculations are carried out within the so-called M- 
scheme where states with M = I are considered. Diago- 
nalizations are done with a parallel shell model program 
we developed a few years ago [47[- Part of the calcula- 
tions are checked with the shell model codes NuShcllX 
Hi and ANTOINE El. All the calculations are done 



is the particle number oper- 



on the computer Kappa at the National Supercomputer 
Center in Linkoping, Sweden. Matrices with dimensions 
up to 10 9 (in the M-scheme) can be diagonalized with 
high efficiency. In Fig. Q]we plotted the M-scheme di- 
mensions for the M = positive-parity states in even- 
even Sn and Te isotopes. The dimensions of the corre- 
sponding 7 W = + states are also given for comparison. 
Only tin isotopes are considered in our fitting procedure 
due to the limitation in computing capability around 
mid-shell when both protons and neutrons are consid- 
ered. 
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FIG. 1. (Color Online) Dimensions of the M 71 = + states in 
even-even Sn (circle) and Te (triangle) isotopes as a function 
of valence neutron numbers, where M and N denote the total 
magnetic quantum number and the number of valence neu- 
trons, respectively. The open symbols connected by dashed 
lines give the dimensions of the corresponding I w — + states. 



The experimental (negative) binding energy of a given 
state i is given by 



^exp =BE cxp (Ar)+Ex(i)) 



(2) 



where BE and Ex denote the binding energy of the nu- 
cleus (with ./V valence nucleons) and the corresponding 
excitation energy of the state relative to the ground state, 
resp ectively. The experimental data are taken from Rcfs. 
[50l IEtI ]. A total number of D = 157 yrast states from 
nuclei 102_132 Sn arc considered. 

We neglect isospin in following discussions for simplic- 
ity since the systems we handle in the present work only 
contain valence neutrons. The calculated total energy of 
the state i can be written as 



E 



cal 



C + Ne + 



N(N - 1) 



V m + (*i\H\*i), (3) 



where "J/ is the calculated shell-model wave function 
of the state i and I is the total angular momentum. 
The constant C denotes the (negative) binding energy 
of the core 100 Sn. The values of £0 and V m depend on 
the way the effective Hamiltonian is constructed. In the 



present work we assume £o 97/2 = V 7 



/2 



and the 
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TABLE I. Numbers of two-body matrix elements (j a jp\V\jsj*t) jt for different sets of J and T. 







J = 


J = 1 


J= 2 


J = 3 


J = 4 


J= 5 


J = 6 


J = 7 


J = 8 


J = 9 


J = 10 


J = 11 


T = 







36 


16 


48 


16 


25 


11 


9 


3 


2 




1 


T = 


1 


15 


6 


46 


18 


34 


13 


16 


6 


4 


1 


1 





other singlc-particlc energies and monopole interactions 
are given as relative values with respect to those of the 
orbital Qg-jn- Thus £o and V m correspond to the real 
energy and monopole interaction of the O57/2 state in 
101 Sn. 

The excitation energy and wave function of a given 
state only depend on the shell model Hamiltonian H. 
One may rewrite the Hamiltonian as H = H m + Hm 
where H m and Hm denote the (diagonal) monopole and 
Multipolc Hamiltonians, respectively. The shell model 
energies can be written as 



E bM = (* 7 |F|* 7 ) 



N a {N p -6 a p) 



(ViIHmI*!), 



where J2 a < N a >= N and 




N a (Np-5ap)\_ N(N-l) 



1 



(4) 



(5) 



We optimize the single-particle energies and monopole 
terms of the realistic effective interaction by minimizing 
the quantity 



(6) 



where the summation runs over all states considered. 
The quality of the fit can be measured by the root-mean 
square deviation as 



X 



D-P' 



(7) 



where P denotes the total number of free terms that are 
to be considered in the fitting. As mentioned before, 
the single-particle ener gy o f the orbital ld§/2 are fixed at 
e(ld 5 / 2 ) = 0.172 MeV 12j. One may rewrite the calcu- 
lated total energy as 



Ef- = £ V k x k + e lds/2 < N ld5/2 > + <H M >, (8) 

k 

where V k denote the unknown singlc-particlc energies and 
monopole interaction terms to be determined. The bind- 
ing energy of the 100 Sn core and the singlc-particlc energy 



eg are also taken as adjustable terms since the uncertain- 
ties in experimental data arc still large [50} . We have 
P = 20 variables in total. 

To minimize the % 2 function we apply a Monte Carlo 
global optimization method which we developed recently 
(denoted as MC in following discussions). It is an it- 
erative approach. The basic idea is as follows: For 
the n-th fitting step we start with an initial set of 
{Vfc(n)} and a new set of variables {V^'(n)} is proposed 
by the Monte Carlo sampling method. We require that 
|V fc '(n) — Vk(n)\ < 8k where 5 k denote the step lengths 
of the sampling. This new set will be accepted as 
{V k {n + 1)} if we have X 2 {Vk( n )) < X 2 {Vk{n)). Another 
consideration is that the proposed new set will also be ac- 
cepted as {Vk(n + 1)} with certain probability g even if 
one has X 2 (^fc( n )) > X 2 (H(^))- This is a key part of the 
global optimization since the x 2 as a function of many 
variables may contain more than one minimum. Other- 
wise the global searching may be trapped in a local one. 
The step is repeated until convergence. The step lengths 
5k and the probability function g can also be adjusted in 
the fitting to getting a faster convergence. 

The advantage of the Monte Carlo global optimization 
method is that no information on the derivatives (i.e., Xk 
in the present study) is required. This is very convenient 
when other observables (e.g., B(E2) values) are included 
in the fitting. 

To get the mean-square deviations x 2 f° r a given set 
of {Vfe(n)} and the large number of succeeding sam- 
plings {VZ(n)} (of the order 10 2 — 10 3 ) one has to rc- 
diagonalizc the corresponding shell-model Hamiltonian 
matrices. Since the shell-model diagonalizations around 
N = 16 arc still time consuming, we further assume that 
the wave functions are stable against the variation of the 
effective Hamiltonian, (4'/(n)|^' / (n)) ~ 1. One has 



p 

Er\n)' - Er\n) « £(V fc '(n) - V k (n))x k , (9) 
fc 



from which the % 2 value for a given sampling can be 
calculated approximately in a straightforward way. The 
wave functions and coefficients x k are re-calculated when 
a new set of variables are accepted. This is known as 
the linear approximation based on which standard fitting 
approaches can be applied [13, HH . 

As a comparison, the singular value decomposition 
(SVD) approach is also employed in the fitting process. 
The SVD approach was used recently in Refs. |53l , l&jj . 
For calculations with the SVD approach, the constants 
C and £0 are taken as free parameters with no rcstric- 
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tion. In the MC approach, we assume that C and Eq can 
only take values within the range —825.5 ~ —8.24 MeV 
and —12 ~ —8 MeV, respectively, by considering the un- 
certainties in experimental data [50j. For the monopole 
interactions we assume \V m \ < 1.5 MeV. These restric- 
tions can be adjusted in the fitting process if necessary. 

The fitting is carried out in three steps. In the first 
step we only consider 131 states in the nuclei 102_112 Sn 
and 12 °- 132 Sn. The nuclei 113 - 114 Sn and 118 - 119 Sn are 
considered in the second step to further fine-tune the 
effective interaction. The three isotopes 115_117 Sn are 
added to the calculation in the last step. Our calculations 
show that convergence is already reached in the second 
step. 

To test the fitting approaches mentioned above, calcu- 
lations are done with 10 sets of random monopole Hamil- 
tonians. They are generated by the Monte Carlo sam- 
pling approach with the restriction 0.172 MeV < e a < 
5 MeV and \V m ., a p\ < 1 MeV. Then these Hamiltoni- 
ans are optimized by fitting to the 131 states mentioned 
above. We found that in all cases one can get conver- 
gence with the MC approach within ten iterations. As a 
typical example, one set of these calculations are plotted 
in Fig. [5] Two types of MC calculations are presented in 
the figure. The solid triangles correspond to the calcu- 
lations with the restrictions on the constants C, £o and 
V m mentioned above. These restrictions are removed in 
calculations marked by the open symbols. This is why 
in the first (n = 1) iteration the x 2 value is smaller in 
the latter case. For n = 1 calculations with the SVD ap- 
proach give the same x 2 value as that of the second MC 
calculation. But the new set of variables {Vk(n = 2)} 
predicted by these approaches are very different. Among 
calculations with the ten random samplings, as in Fig. [2] 
the SVD approach diverges in most cases. The reason for 
the divergence may be that the step lengths (the differ- 
ence between {I4(n + 1)} and {T4(n)}) predicted by the 
SVD approach is too large that one can not apply the lin- 
ear approximation. In the MC approach we restrict the 
step length be 6 < 1 MeV. This is a rational restriction 
by taking into account the fact the T = 1 monopole inter- 
actions are mostely small and close to zero in shell-model 
calculations in light and medium- mass nuclei [35l l38j . 

Fig. [2] suggests that our restrictions on the C, £o and 
V m are reasonable and have no influence on the final re- 
sults. It should also be mentioned that the monopole 
Hamiltonians that arc fitted starting from these random 
samplings are very similar to each other. 

In Fig. |3] the same calculations are done starting from 
the realistic CD-Bonn potential. The 131 yrast states are 
considered in the fitting. The initial root-mean-squarc 
deviation is around 750 keV. This is reduced to about 
110 keV after 10 iterations. The x 2 given by the SVD 
approach is smaller than the MC method for the same 
reason as before. But the variables C, £o an d V m pre- 
dicted by the SVD approach soon fall into the expected 
range after a few iterations. The x 2 value given by the 
SVD calculation starts to fluctuate around n = 10. This 
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FIG. 2. (Color Online) The convergence of the mean- 
square deviations \ 2 as a function of iteration number within 
the singular value decomposition (SVD) and Monte Carlo 
global optimization (MC) approaches starting from a random 
monopole Hamiltonian. See text for the difference between 
the red-solid and green-open triangles. 



is avoided in the MC calculations by gradually decreas- 
ing the maximal range 8. In the calculation labeled by 
the open circle presented in the figure, the MC approach 
is applied instead of SVD. 

In Ref. [55| the maximal deviation r = max|S| al — 
E^ xp | is minimized in their fitting to experimental bind- 
ing energies within the Skyrme-Hartree-Fock approach. 
In the insert of Fig. [3] we plot the maximal deviation r. 
It also gradually decreases as a function of the iteration 
even though the criterion of our calculation is to find the 
minimum of x 2 ■ 
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FIG. 3. (Color Online) The convergence of the mean-square 
deviations x 2 as a function of iteration number for fittings 
that are started from the realistic CD-Bonn interaction. The 
insert shows the maximal deviation r = max|£^ al — -E^ xp | in 
each step. 
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TABLE II. The final mean square deviation \ 2 ( in MeV 2 ) 
given by the SVD and MC fittings and the predicted values of 
the variables C, so and V m and their uncertainties (in MeV). 

X 2 C £0 Vm 

MC 2.697 -825.5 ±0.288 -10.666 ± 0.0410 0.167 ± 0.00249 
SVD 2.686 -825.5 ±0.288 -10.669 ± 0.0410 0.172 ± 0.00249 



III. RESULTS AND DISCUSSIONS 

The fitting approach described above has been success- 
fully applied in deriving the effective interactions for sev- 
eral nuclear regions. In particular, we constructed a new 
interaction to describe the structure of the TV = 83 — 85 
isotones by assuming the nucleus 146 Gd as the core [5l| . 
which is not as good shell closure as 100 Sn. As a result, 
the standard SVD approach failed to get a converging 
result. The new interaction for light tin isotopes derived 
in the present work is briefly discussed below. It has al- 
ready been used in a recent experimental study on the 
level structure and electromagnetic transition properties 
of the odd- A nucleus 109 Te 

The monopolc Hamiltonians we derived in Fig. [3] arc 
slightly refined by including the nuclei 113_119 Sn into the 
fitting. As mentioned before, we include a total num- 
ber of 157 states in the fitting among which there are 
31 binding energies and 126 excitation energies. As seen 
in Table ITI1 after around 15 iterations both calculations 
give a mean-square deviation \ 2 ~ 2.7 MeV. It means 
that these states can be reproduced within an average 
deviation of about 130 keV. In Table llTl we also give the 
values of the variables C, eq and V m predicted by the MC 
and SVD calculations. The uncertainties within these 
variables are also analyzed with the help of the SVD ap- 
proach. The C and e values predicted by the fitting are 
in reasonable agreement with the binding energies given 
in Ref. g^. 

The largest uncertainties of the optimized monopolc 
Hamiltonian are related to the single-particle energies. 
The values predicted by the MC approach are £id 3/2 = 
5.013 ± 3.10, e 2si/2 = 0.369 ± 2.63 and e 0hll/2 = 3.249 ± 
0.83 MeV. The single-particle energies given by the SVD 
approach are £id 3/2 = 5.182 ± 3.16, £2s 1/2 = 0.409 ± 
2.73 and e h 11/2 = 3.236 ±0.84 MeV. More experimental 
efforts are desired in order to get a better constraint on 
these single-particle energies. 

In Table |HT1 we compare the optimized monopole terms 
from the SVD and MC approaches and those of the re- 
alistic CD-Bonn interaction. The list of the two-body 
matrix elements and part of the calculated results on tin 
isotopes with the optimized interaction can be found in 
Ref. IH. 

The monopole Hamiltonians optimized with the SVD 
and MC approaches are similar to each other. To illus- 
trate this point, in Fig. 0]we plotted the deviation from 
experimental data for calculations with the two effective 



TABLE III. The strengths of the T = 1 monopole interactions 
Vj^-jj (in MeV) of the original CD-Bonn potential and the 
interactions optimized by the SVD and Monte-Carlo fitting 
approach. 

ji j2 CD-Bonn SVD MC 
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FIG. 4. (Color Online) Differences between experimental and 
calculated binding energies, £^ xp — Ef al , as a function of va- 
lence neutron number. 



Hamiltonians. The difference between the calculations is 
practically negligible. The largest deviation from exper- 
iments is seen at N = 15 where the energy 5/2+ state 
is over-estimated by around 600 keV. In the following 
only calculations done with the MC optimized effective 
Hamiltonian arc presented for simplicity. 

The calculated shell model energies, E SM =< H >, 
of the selected states in tin isotopes are plotted in Fig. 
[5] The contributions from the monopole Hamiltonian are 
also presented for comparison. From the figure one can 
see that the contribution from the multipole Hamiltonian 
reaches its maximum around the mid-shell. 

The excitation energies of the 126 excited states can 
be reproduced within an average deviation of about 180 
keV. The largest difference is seen in the 5/2+ state in 
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FIG. 5. (Color Online) The calculated shell model energies 
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FIG. 6. (Color Online) Calculated excitation energies of the 
low-lying even-spin states in 102 ~ 130 Sn. 



115 Sn mentioned above, where the experimental datum is 
under-estimated by about 1 MeV. Extensive experimen- 
tal efforts are made recently to explore higher-seni ority 
states built on the 10 + states in Sn isotopes [13, EH, l59| . 
In Fig. |5] we plot the calculated excitation energies of 
the low-lying even-spin states up to I 11 = 12+ in the nu- 



cle 



: 102-130 



Sn. The experimental data are n ot p lotted for 



simplicity. They can be found in Refs. [171 Il8l. l51l. l59j. 

In Fig. [7] we present the calculated occupancies of the 
single-particle orbitals in the ground state wave functions 
of even tin isotopes. The structure of the low-lying states 
in light tin isotopes are dominated by configuration mix- 
ing between the orbitals 0g 7 /2 and 1^5/2- The two or- 
bitals are half-filled around N = 8 ( 108 Sn). The orbital 
0^11/2 is calculated to be half-filled around N = 22 and 
24. This is in agreement with the speculation in Ref. 
[60| by considering the E2 decay properties of the 10 + 
isomers in these nuclei. 

The calculated excitation energy of the 11/2 - , l/2 + , 
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FIG. 7. (Color Online) Shell model occupancies, < Nj > 
/(2j + 1), of the three higher- j shells in the ground states of 
even tin isotopes. 



5/2 + and 7/2 + one-neutron-hole states relative to the 
3/2+ state arc -0.022, 0.475, 1.720 and 2.521 MeV, re- 
spectively. These are in fair agreement with the experi- 
mental data I5l], roll . 



IV. SPIN INVERSION IN llu SN 



The 
state i 



33 
12 



spins of the ground state and first excited 
1 103 Sn are / = 5/2 and 7/2, respectively 
which are reversed with respect to those in 101 Sn 
Through seniority model analyses with a pair- 
ing Hamiltonian, Ref. 12| suggested that the inver- 
sion is dominated by orbital-dependent pairing corre- 
lations, namely the strength of the pairing matrix el- 
ements 



W,2 



(ld 2 5/2 \V\ldi /2 ) J= 



V\0g 7/2 



)j=o is much larger than that of 
This produces strong additional 
binding for the = 5/2 + state in 103 Sn, which even- 
tually becomes the ground state. The effect of other in- 
teraction terms on the spin inversion was not considered 
in Ref. [13]. 

In Fig. [5] we analyze the contribution from different 
components of the effective Hamiltonian to the spin in- 
version between 101 Sn and 103 Sn. For that purpose the 
energies of the first = 5/2+ and 7/2 + states are cal- 
culated with a limited Hamiltonian H' containing the 
single-particle terms and two-body matrix elements with 
J < 4ax only. J max denotes the maximal spin value 
of the two-body matrix elements to be considered. Two 
different types of calculations are presented in the figure. 
The solid symbols correspond to the results calculated 
by diagonalizing the Hamiltonian H' . The expectation 
values of such a Hamiltonian with respect to the corre- 
sponding wave functions \^i) of the full Hamiltonian H, 
(^fi\H'\^fi), are plotted as open symbols. It should be 
mentioned that calculations with the original CD-Bonn 
interaction and other realistic nuclcon-nuclcon potentials 
give a similar pattern. 
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FIG. 8. (Color Online) The calculated shell-model energies of 
the first 5/2 + and 7/2 + states in 103 Sn with a Hamiltonian 
H' containing the single-particle terms and two-body matrix 
elements with J < J max - The solid symbols denote the eigen- 
values of H' derived by diagonalizing the Hamiltonian ma- 
trix. The open symbols denote the expectation values of H', 
{^i\H'\^/i}, with respect to the corresponding wave functions 
of the full Hamiltonian H , where all matrix elements are 
taken into account. 



It is thus seen from Fig. [8] that both calculations give 
similar results concerning the order of the 5/2+ and 7/2 + 
states. Calculations with the pairing matrix elements 
only (i.e., J ma x = 0) show that the pairing terms, in par- 
ticular the element (0gy 2 \V\0gy 2 ) j=o mentioned above, 
can significantly reduce the gap between the two states 
but were not strong enough to induce the inversion. A 
sudden switch is seen when the J = 6 two-body matrix 
elements are considered. It can be seen from Fig. [5] that 
in both calculations the exact results are also approached 
by including terms with J < 6 only. This is expected 
since the low-lying states of light tin isotopes mainly oc- 
cupy the nearly degenerate orbitals 0<?7/ 2 and ld 5 / 2 for 
which the maximal spin is J — 6. Thus the contribution 
from interactions with higher spin values are marginal. 
Among T = 1 matrix elements the maximal spin one can 
have is J = 10. It corresponds to the coupling of two 



nucleons in the orbital 0hn/ 2 . Calculations in the re- 
stricted 0(77/21(^5/2 model space give a result similar to 
Fig. 1 

Among the J = 6 two-body matrix elements the ones 
that contribute most to the spin inversion are the repul- 
sive matrix element (0gy 2 \V\0gy 2 ) j=q and the strongly 

attractive one {QgT /2^-d^/ 2 \V\Qg-j /2^d^/ 2 ) j=%- This can 
be understood by considering the structure of the wave 
functions of the two states. In the 5/2 + ground state in 
103 Sn, the leading component is \(0g 7 / 2 ) j=o^d 5 / 2 ) i E3- 
Its overlap with the total wave function is calculated to 
be \((0gy 2 )j= ld 5 / 2 \^}i\ = 0.86. One can also construct 
a three-body state starting from the pair \0gr/ 2 ld^/ 2 )j = e. 
The overlap between the state thus constructed and the 
total wave function is |(0g7/ 2 ±G? 5 / 2 )j = 60g7/ 2 | , I')/ = 0.75. 
The (0g 7 / 2 ld 5 / 2 \V\0g 7 / 2 ld 5 / 2 ) J=6 term induces a sig- 



nificant additional binding for the J 71 — 5/2 + state 
in 103 Sn, as can be seen from Fig. [5] It should 
be mentioned that states generated by the two cou- 
plings |(037/ 2 )j = olrf 5 /2} and (0g 7 / 2 ld 5/2 ) J=6 0g 7 / 2 ) i arc 
not perpendicular to each other. Their overlap is quite 
large, ((0g2 /2 ) j=0 ld 5/2 |(0 ff 7/ 2 ld 5 / 2 )j =6 037/2}j = 0.74. 
This can be evaluated analytically [62l |. 

The overlaps of the total wave function of the first 7/2 + 
state 103 Sn with its leading components are calculated to 
be K(ld» /2 )j= 5r/2 |*)/|=0.65 J |<(0fl? /a )./05r /2 |*)/| = 
0.62 and | <(0sr 7/2 ld 5/2 ) , 7=6 ld 5/2 / 1=0.57. Fromashell- 
model point of view, the couplings \(0gy 2 )j=o0g 7 / 2 )i 
and | (0^7 /2 ) J=6 0<?7 /2 ) I generate exactly the same three- 
particle state. All interaction terms {0gy 2 \V\0g 2 /2 )j 
contribute to the total energy of the state [631. It ma y 
be interesting to mention that the effect of the maxi- 
mally aligned pair in single- j systems was discussed in 
Ref. 111,0. 




FIG. 9. (Color Online) Solid symbol: The overlaps between 
the wave functions \9i) of the full Hamiltonian H and those of 
the pairing Hamiltonian with J max = for the first 5/2 + and 
7/2 + states in light odd-A Sn isotopes. Open symbol: Same 
as above but only the non-diagonal pairing matrix elements 
are considered. 



The influence of the pairing interaction on the struc- 
ture of tin isotopes was considered in a variety of ap- 
proaches (see, e.g., [27], HU). In Fig. [S]we calculated the 
overlaps between the full shell-model wave functions and 
those calculated from the pairing Hamiltonian including 
the single-particle terms and the J = pairing matrix 
elements. It is thus seen that the non-diagonal pairing 
matrix elements (j 2 \V\j' 2 ) j = o play an essential role in 
inducing the configuration mixing in the first 5/2 + and 
7/2+ states of odd-^4 tin isotopes. The overlap gradually 
decreases when the number of valence neutrons increases. 
This is consistent with the generalized seniority model 
calculations in Ref. [26j], namely the overlap between 
the full wave functions and seniority-truncated state also 
decreases with increasing neutron number. 



8 



V. SUMMARY 

The structure properties of light tin isotopes are calcu- 
lated with a global optimized effective interaction. The 
unknown single-particle energies of the orbitals ld 3 / 2 , 
2si/2 and O/111/2 and the monopole interactions are re- 
fined by fitting to experimental binding energies. A total 
number of 157 states in 102_132 Sn are considered in the 
fitting. The binding energies of these states can be re- 
produced within an average deviation of about 130 keV. 
The largest deviation is around 600 keV which is seen 
in nuclei 115 Sn. With the effective Hamiltonian thus de- 
rived we calculate the contributions of the monopole and 
multipolc Hamiltonian to the total shell-model energies. 
The excitation energies of the low-lying even spin states 
in even Sn isotopes are presented. We also evaluated the 
shell model occupancies in the ground states of these nu- 
clei. Detailed systematic calculations on the spectra and 
decay properties of tin isotopes as well as the list of the 
two-body matrix elements will be presented in Ref . [58| . 



We analyze the origin of the spin inversion between the 
7/2 + and 5/2+ states in 103 Sn and heavier odd tin iso- 
topes in order to explore the possible influence of different 
interaction channels. We thus find that both the J = 
pairing and the maximally aligned J = 6 two-body ma- 
trix elements produce strong additional binding for the 
5/2 + states. The non-diagonal pairing matrix elements 
play an essential role in inducing the mixing of different 
configurations in the wave functions of these states. 
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